Current without bias and diode effect in shuttling transport of nanoshafts 
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A row of parallely ordered and coupled molecular nanoshafts is shown to develop a shuttling 
transport of charges at finite temperature. The appearance of a current without applying an external 
bias voltage is reported as well as a natural diode effect allowing unidirectional charge transport 
along one field direction while blocking the opposite direction. The zero-bias voltage current appears 
above a threshold of initial thermal and/or dislocation energy. 
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Organic field effect transistors (OFETs) based on dif- 
ferent polymers [l|, [2|, |3[ attract an increased interest due 
to numerous potential applications as flexible and low 
cost storage and microelectronic devices. For achieving 
excellent electric properties such as high charge carrier 
mobilities and low resistive losses, which are required 
for technologically attractive device applications, a high 
structural ordering of the semiconductor molecular mate- 
rial is necessary [4J, [5| . Oligothiophenes and their deriva- 
tives can be regarded as one of the most promising sys- 
tems for building such self-organized structures across 
multiple length scales due to the variety of intra- and in- 
ter molecular interactions, which originates from the po- 
larizability of the sulfur electrons and the aromatic tt- 
electron system [6|, [D, |8[ . Just recently an OFET struc- 
ture has been built from ultra-thin self-assembly films 
made up from oligothiophenes, which are arranged in a 
high-order lamellar stacking perpendicular to the sub- 
strate surface [3|. 

In general, the charge transport is largest in the di- 
rection perpendicular to the plane of the thiophene rings 
|Sl, LLQ[. This finding addresses already the basic differ- 
ence between electronic transport in organic conductors 
and in classical semiconductors. Though the band gap 
between the lowest unoccupied and the highest occupied 
level of about 3 eV in quarterthiophene suggests an anal- 
ogy to the conduction and valence band in semiconduc- 
tors, there are crucial differences. While in conventional 
semiconductors the transport is due to delocalized states 
and limited by the scattering, the transport in molecules 
is due to localized states dominated by hopping [111 ]. 

In this paper we suggest a new mechanism of charge 
transport which is possible for flexible molecular tubes, 
shafts or any elastically deformable assembly which could 
provide an alternative explanation of transport proper- 
ties in nanoshaft-based OFETs. The performance of 
OFETs based on oligothiophene |3[ shows characteris- 
tic features, e.g. the current starts at a certain threshold 
of gate voltage and reaches a saturation value for certain 
drain voltages. We will show in the following that for a 
regular arrangement of long elastic molecules a shuttling 



transport can be established, which can model these fea- 
tures. In particular, we will describe a shuttling effect 
of nanoshafts leading to a diode behavior as well as to 
a directed current without external bias voltage. This 
resembles the ratchet effect described as Brownian mo- 



tors h^,[l3|,[l4| also realizable with oscillating laser fields 
15l llq . The thermal noise plays an important role for 



activating such motors to overcome a certain barrier [17[ . 
In the effect presented here the activating threshold is 
given by an initial kinetic or potential energy necessary 
for the first nanoshaft to reach the contacts. 

Shuttling in single-dot devices had been already in 
the center of interest [18|, [l9[. Besides promising ap- 
plications for nano-electromechanical devices the cou- 
pling of quantized transport with micromechanical can- 
tilevers [20f] bridges the classical mechanical and the 
quantum physics in an exciting way. Our model, ana- 
lytically solvable, demonstrates that a similar shuttling 
mode as known from single-dot transport by mechani- 
cal cantilevers is achievable by a row of elastically cou- 
pled molecular tubes. The model is based on the coupled 
spring chain. Despite their simplicity, such models sur- 
prise with the ability to describe unknown phenomena of 
classical quantization phenomena in velocities 21]. The 
model of freely rotating linear chains has been investi- 
gated also in view of quantizatio n 12 2| and has been ap- 
plied to entanglement problems [23[. We note that the 
shuttling transport described here should not be mixed 
up with the so called 'chain shuttling' in polymerization 
reactions, where a growing polymer chain is transferred 
between different metal catalysts [24|, [25| . 

The elasticity of polymers with respect to the bending 
and rigidity has been studied by molecular dynamic sim- 
ulations [26]. The model of hard-sphere chains is used 
to describe the stiffness and diffusion in polymers [271 ] . 
These elastic properties play an important role in the 
mesoscale modeling of carbon nanotubes [28[ . Such car- 
bon nanotubes on glassy carbon electrodes can enhance 
the current response by a factor of 1000 as reported in 
[29[. Therefore it is of high actual interest to reinves- 
tigate the transport in nanotubes. We will concentrate 
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FIG. 1: Snapshots of time evolution of the chain of shuttling nanoshafts (left) with positive charge (green) and negative charges 
(red) or no charge (grey). The kinetic energy, currents on both sides and the recombination rate are plotted from left to right. 



in this paper on charge transport perpendicular to the 
direction of nanoshafts. 

Specifically, we consider N nanoshafts of length L par- 
allely oriented in the ^/-direction. The shafts are standing 
perpendicular on the substrate, and are periodically ar- 
ranged along the x-direction at a certain distance from 
each other (see Fig. [1]). We expect them to couple elas- 
tically by their top ends. Then the shafts bend in x 
direction when applying the tangential force F with the 
displacement 
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where 6 is the Young modulus and / = J z 2 dA is the 
area moment of inertia being / = Trr^/A for circular 
cross sections of tubes with radius ro or / = ab 3 / 12 for 
nanoshafts with rectangular cross section of sides a and 
b. The maximal displacement on the top ends z = L 
is x = x(l) = FL 3 /3£I from which the spring con- 
stant k = F/x = 361 /L 3 is given provided we know 
the Young modulus. In other words we have a coupled 
linear chain of top ends obeying the differential equations 
■^Xi = oo 2 (xi-i — 2x{ + Xi+i) with uj 2 = k/m given by 
the mass of the nanoshafts. The nanoshafts are located 
in-between two solid contacts such that the coupled equa- 
tions are bounded which can be formally expressed by 
xo = £jv+i = 0- Additionally each nanoshaft can carry 
a charge qi which can be exposed to a time-dependent 
external field E(t). The total coupled equation system 
for the time-dependence of the nanoshaft top ends reads 



d 2 Xj(t) 
dt 2 



■ UJ 2 \Xj- 



l (i)-2x i (t)+x i+1 (t)]+a i (t) (2) 



with uj 2 = 361 jvnL 3 and ai(t) = qiE(t)/m. For simplic- 
ity we will use the time in units of 1/u further on. For 
quarterthiophene we have a typical value of k = lN/m 



and a mass of m = 5.49 x 10 _22 g which leads to the scale 
uo = 1.35/ps. The length will be in units of /q = O.lnm 
as typical value. The field is given in units of mcj 2 lo/q = 
6.2 x 10 8 V/m again given explicitly for quarterthiophene 
and the charge current in units of quo — 21.6mA where 
we assume one single charge q = e. The typical energy 
units are mu 2 ll/2 = 31.2meV= 362K. 

The analytical solution of ([2|) obtained by the standard 
eigenvalue method is given in terms of the normalized 
orthogonal system (j) nv = y /f 2/(N + 1) sin [nu7r/(N + 1)], 



N 



m = E' 



n=l 



(c n cos uj n t + d n sin uo n t) 



z 

I 



dt 



, sin cu n (t — t') 



N 



J2 a m(t') 



(3) 



with the eigenfrequencies uj 2 = 2(1 — cos [nir / ( (N + 1)]) 
and the initial condition determining c n = ^2 U 4>nuX u (0) 
and d n = Y, u <t>nvX v {ti) / u n . 

Now we proceed by considering the nanoshafts between 
two oppositely charged plates, g S ides? modeling the con- 
tacts. Each nanoshaft can carry a negative charge qi < 
describing the charge in the lowest unoccupied molecular 
orbital or a positive charge qi > for the transport of 
holes in the highest occupied molecular orbital or none. 
Each time when two top ends of the nanoshafts touch 
each other the charge is moved if one of both tubes had 
no charge. In the case of opposite charges they annihi- 
late and it is counted as recombination which gives rise 
to light emission. If one nanoshaft is touching the con- 
tact on one side it contributes to the current on this side 
as Asides — Qi • Each time such an event happens the time 
evolution according to the analytical solution (|3]) restarts 
with the new initial conditions and new charge distribu- 
tion. In this way we use the speed of analytical solution 



together with the nonlinear process of recharging. 

For exploratory reasons we restrict ourselves to 11 
nanoshafts of length 20/o and have arranged them at a 
distance of 3/o- These are typical parameters of quar- 
terthiophene used in recent experiments [3[. First we 
charge all shafts equivalently and start the process of 
shuttling by bending the most left tube to the left con- 
tact, allowing to transfer the first charge to the left lead. 
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FIG. 2: The time evolution for the recombination rate 
(above), the total kinetic energy (middle) and the total cur- 
rent (below) for 3 different applied fields. The time unites are 
in cu 2 = 3EI/mL 3 and the length units are /o- 

In figure [T] we have plotted different snapshots of the 
time evolution in units of 1/uj. We see that the system 
starts to shuttle and to transport positive charges from 
the left to the right and negative charges from the right to 
the left side. This leads to steps in the current counted as 
charges delivered on the corresponding side per elapsed 
time. The recombination rate is increasing with time; 
complete movies can be found in [301 ]. 

The snapshots in figured] are actually taken from a run 
without applying any external field. The gates serves 
here merely as a reservoir of charges, not as capacitor 
plates. The astonishing observation is now that even for 
such an unbiased case a finite total current is developing. 
This fact is demonstrated in figure [2] where the time evo- 
lution of the recombination rate, the kinetic energy, and 
the total current which is the sum of right and left cur- 
rents, are plotted for different external fields. The total 
current and the recombination rate reach saturation with 
increasing time while the kinetic energy remains constant 
on average. If one applies an electric field of E = 0.2 op- 
positely to the direction of zero bias voltage current we 
see that the kinetic energy is decreasing with time and 
the shuttling current as well as the recombination rate 
are reduced. In fact, it reaches a state where the left and 
right half of the tubes shuttle symmetrically against each 
other such that no net current is delivered. If the field 
is applied in the opposite direction E = —0.2, the sys- 
tem is accelerated and delivering more and more charges, 



limited finally only by the length of the tubes. In other 
words we observe a pure diode effect just due to shut- 
tling of coupled chains of nanoshafts. The total current 
reached after long times for different applied fields is de- 
picted in figure [3] and illustrates this diode effect as a 
quite systematic one. 

One may argue that this diode effect as well as the 
observation of a current without external bias is due to 
the initial charging of the nanoshafts. In the figure [3] 
we have plotted the finally reached total current versus 
the applied field for two cases, the initially charged case 
considered so far, and the case where initially all shifts 
are uncharged. We see that identical final currents are 
obtained without external bias for both cases and the 
diode effect is present, as well. 
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FIG. 3: The total current versus applied electric field for two 
cases, initially uncharged tubes and initially charged tubes. 



Now we want to return to the puzzling observation that 
a charge transport occurs even without external bias. To 
understand the conditions under which such an effect oc- 
curs we have to clarify the dependencies on the initial 
kinetic and potential energy and the geometric config- 
urations. The length L of the nanoshafts determines 
only the time scale uj and scales out here. Therefore 
the distance between the nanoshafts Ax in the row is 
the only geometric parameter left. Further, the initial 
kinetic and potential energy we determine by the elonga- 
tion xq and speed ^o of the leftmost nanoshaft while all 
other nanoshafts are initially not displaced and at rest. 
In figure 0] we plot the dependence of the reached cur- 
rent on the initial displacement and the initial velocity. 
Using dimensionless units for the velocity v* = v/uoIq 
and the displacement xo* = xq/Iq the initial energy is 
^tot = mu 2 lQ(vo * 2 +£o* 2 )/2 where the length unit is 
/o- We see in the right plot of figure [4] that the cur- 
rent is increasing quadratically with the initial velocity, 
i.e. the current is proportional to the initially deposited 
kinetic energy. This initial kinetic energy can be equiva- 
lently realized by thermal motion and characterized by a 



temperature. In the simulation described so far we bend 
the leftmost nanoshaft sufficiently close to the electrode 
such that the shuttling can happen even without initial 
kinetic energy. If we bend it less we obtain a threshold 
in kinetic energy analogously to the thermal threshold in 
Brownian motors [17]. The dependence of the current on 
the initial position is seen in the left plot of figure [U We 
have chosen an initial velocity below the threshold which 
shows that a threshold in the position or potential energy 
has to be overcome in order to create the shuttling cur- 
rent. We note that the shuttling transport develops only 
if we start with asymmetric initial conditions. Otherwise 
counter-oscillations of the left and right half of shafts 
block any transport. Therefore the symmetry breaking 
is due to initial conditions and this model illustrates an 
electronic ratchet effect with thresholds in kinetic or po- 
tential energy corresponding to a high system tempera- 
ture or initially elongation being mcu 2 Ax 2 /2 = 281meV 
for the here used example of quarterthiophene. 
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FIG. 4: The total current without voltage bias versus initial 
displacement (left) and versus initial velocity (right) of the 
most left nanoshaft. 

In summary we have found that a chain of perpendic- 
ularly arranged coupled and chargeable nanoshafts show 
a shuttling transport of charges. As a surprising effect 
it turns out that a finite current is established already 
without external bias only due to the initial asymmet- 
ric deformation of the nanoshafts. The effect observed 
here is reminiscent of the ratchet effect. Even in the ab- 
sence of a net macroscopic force or noise a current can be 
generated. This resembles a lot the prominent motor pro- 
teins found inside cells. There, proteins such as Mysosin 
and Kynesin use chemical energy which is gained from 
hydrolysis of ATP into ADP to move along asymmetric 
pathways transporting vesicles inside cells, contracting 
muscles and being important in the process of cell divi- 
sion 
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